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Sucasary 

In  a  aeries  of  reoent  papers  the  prlnolple  of  invariant 
imbedding,  which  xrepresents  an  extension  of  the  invariance 
principles  of  Ambarsumian  and  Chandrasekhar,  has  been 
ecqployed  in  the  study  of  a  variety  of  physical  processes 
including  radiative  transfer,  neutron  transport,  random  walk 
and  scattering,  and  wave  propagation. 

use  of  this  principle  leads  to  the  formulation  of 
various  fimctional  eqviations  describing  the  process  under 
consideration.  Redheffer  has  used  a  similar  approach  in  the  * 
treatment  of  several  problems  of  electromagnetic  theory. 

The  aim  of  the  present  paper  is  to  extend  previous  results 
and  techniques  to  cover  cases  involving  plane  wave  propagation 
in  stochastic  media. 

In  general  terms,  the  approach  Involves  first  the 
derivation  of  stochastic  functional  equations  for  reflection 
and  transmission  coefficients,  followed  by  the  taking  of 
esqpected  values  of  appropriate  functions  of  the  random  state 
variables.  This  makes  possible  the  determination  of  their 
characteristic  functions  and  distribution  functions,  by  swans 
of  still  other  functional  equations,  or  by  computational 
schemes  of  the  Monte  Carlo  type. 

We  discuss  the  particular  example  In  which  a  plane  wave 
is  incident  on  a  stratified  slab  which  is  characterized  by 
stochastic  wave  numbers  in  each  stratum.  !nie  distribution 
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fwiatlous  for  the  ai^litude  of  the  random  reflected  and  trans¬ 
mitted  vfaves  are  then  determined  as  functions  of  the  thickness 
of  the  slab.  The  effects  of  omltiple  scattering  are  taken 
into  account,, 

The  emphasis  thx*oughout  is  upon  the  general  applicability 
of  these  functional  equation  procedures  to  problems  of  propa¬ 
gation  in  stochastic  media,  using  methods  quite  distinct  from 
the  classical  techniques. 
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INVARIANT  IMBEDDING  AND  WAVE  PROPAGATION 
IN  STOCHASTIC  MEDIA 

Hlohaxd  BslXman  and  JRobtrt  K&laba 
Tha  HAND  Corporation^  Santa  Monioa^  California 

1 ♦  Introduction 

In  this  paper,  w«  vrish  to  describe  an  £^plloatlon  of  the 
theory  of  invariant  labedding,  [4  «  llj  ,  to  the  propagation 
of  eleotronagnetio  waves  throui^  an  inhonogeneous  sedim, 
especially  where  the  inhoinogeneity  is  stochastic  in  origin, 
problems  of  this  sort  arise  in  a  variety  of  fields  including 
radio  wave  propagation  and  acoustics. 

There  have  been  a  large  nualber  of  papers  devoted  to 
this  topic,  most  using  the  routine  perturbation  approach,  and 
a  few  the  concept  of  a  stratified  medium;  cf.  Redheffer,  > 
Luneberg,  [l9]  #  Eremaer,  [l3]  #  Tatarskii,  [24 j ,  Gronwall,  [^^j, 
et  al  t  Topics  of  current  interest  in  radio  wave  propagation 
are  discussed  in  [l2,l4,25#27j . 

Our  aim  is  first  of  all  to  illustrate  the  applicability 
of  functional  equations  and  principles  of  invariance  to  the 
study  of  various  types  of  wave  propagation.  The  basic  ideas 
were  sketched  in  our  note,  [7]/  bjoA  fuller  mathematical 
details  will  be  given  in  papers  shortly  to  appear.  Secondly, 
we  wish  to  point  out  the  isq^ortanoe  of  functional  equations  ' 
in  stochastic  variables,  prior  to  the  appearance  of  any 

- s - 

In  his  book  on  Radiative  Transfer,  Chandrasekhar  refers 
to  some  earlier  uses  of  invariance  principles  by  Stokes  and 
Rayleigh. 


Qxpeotad  values.  'Hiftsa  relations  can  be  used,  as  discussed 
below,  to  furnish  nuioerlc&l  solutions  ailong  Ncnte  Carlo  lines. 
Finally,  we  wish  to  en^jhasize  that  the  nonlinear  aspects  of 
the  recurrence  relations  that  are  derived  medce  it  ia?)ortant 
not  to  use  expected  values,  but  rather  to  examine  the  actual 
probability  distribution  of  the  random  variables  that  appear, 

k*  Physical  Prooesa 

tfe  wish  to  consider  the  problem  of  determining  the 
oharaoteristios  of  waves  reflected  from  a  randomly  Inhomo^ 
geneous  medium  and  the  properties  of  waves  transmitted 
through  such  a  nediiua.  In  particular,  we  limit  ourselves  to 
the  case  in  which  a  plane  scalar  wave  is  normally  incident 
on  a  slab  bounded  by  parallel  planes,  the  wave  number  or  index 
of  refraction  of  which  is  a  random  function  of  distance  from 
an  interface.  We  assume  that  the  media  in  which  the  slab  is 
ccmtained  have  constant  wave  numbers  and,  as  usual,  that  time 
variations  are  sisqple  harmonic.  Ultimately  we  are  concerned 
with  random  solutions  of  the  reduced  wave  equation  in  one 
dimension 

(1)  w”(x)  +  k^(x)w(x)  o  0, 

where  the  local  wave  number,  k(x),  is  a  random  function  of 
X.  The  index  of  refraction,  n(x),  is  ©xpreesable  as 
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Since  it  is  the  method  of  attack  which  we  wish  to 

emphasize,  before  taking  up  the  stochastic  case,^  we  solve  the 

reflection  and  transmission  problems  for  the  case  in  which 

k(x)  is  deterministic,  making  use  of  the  notion  of  Invariant 

imbedding  and  a  wave- localization  principle  which  we  first 

» 

gave  in  .  Then  we  take  up  the  problem  of  reflections  from, 
and  transmissions  throxigh,  a  medim  which  has  z^dom  stx^ati- 
fications,  each  stratum  being  bounded  by  planes  perpondloular 
to  the  direction  of  propagation  of  the  incident  wave  and  each 
having  a  random  wave  number.  For  oonvenience,  we  assusw  that 
these  random  wave  numbers  are  Independent,  but  no  assus^tions 
are  made  concerning  "staallness"  of  departures  of  the  wave 
numbers  from  their  average  values. 

Some  generalizations  and  indications  of  future  work  are 
provided  in  the  discussion  of  ^6. 

3.  Deterministic  Wave  Propagation 

To  illustrate  some  of  the  ideas  which  we  shall  en^loy  in 
the  treatment  of  the  stochastic  problem,  we  first  discuss  the 
problem  of  the  reflection  of  a  plane  wave  nonaally  incident 
on  an  inhomogeneous  slab,  a  problem  of  some  difficulty  and 
l^ortanoe  in  itself,  ^19]  .  Consider  a  plane  wave, 
exp  i(kQX  — at),  arriving  from  the  homogeneous  space,  x  <  0, 
and  approaching  the  inhomogeneoi^  slab,  0  <  x  <  b.  !0ie 
space  X  >  b  is  assumed  to  be  homogeneous. 
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Pig.  1 

A  Plane  Wave  Inoident  on  an  InhoiaogeneouB  Slab 

The  mediiia  x  <  0  Is  oharaote^sed  by  the  wave  nuiober 
and  the  medlm  x  >  b  by  the  tfave  number  k, .  The  Inhoao- 

Cm* 

geneoxis  slab,  0<x<b,  is  oharacterlzed  by  the  wave 

nujnber  k(x),  which  Is  assumed  to  be  sectlonally  smooth. 

Our  objective  is  the  determination  of  the  reflected  wave, 
-ik^x 

Ae  ^  ,  where  here,  and  in  what  follows,  we  shall  suppress 
the  time  factor  e”^^. 

We  first  imbed  this  problem  within  a  class  of  problems 

in  which  the  incident  wave  md  the  inhomogeneous  slab  are 

suitably  altered,  f^re  precisely,  we  consider  the  problem 

of  determining  the  wave  reflected  In  the  half-space  x  ^  z, 

where  0  ^  z  b,  as  a  result  of  having  the  wave 

exp[ik;(z  —  0)(x  —  z)j  incident  on  the  inhomogeneous  slab 

z  ^  X  ^  b.  As  usual,  k(z  —  0)  denotes  the  limit  of  k(x) 

as  X  approaches  z  from  below  and  k(z  f  0)  is  the  limit 
* 

as  X  approaches  z  from  above.  As  before,  we  take  the 
space  b  <  X  to  be  homogeneous  with  che  wave  number  k^^^. 
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k(z  —  0) 

-  - . -1  - . . . . 

lt(x) 

h 

1 

0  1 

t  1 

1  1  r  ^ 

Plg.  2 

The  O«oi&«tzy  of  tho  Class  of  Slab  Problaias 


For  tho  laodifltd  probltm  th«  spaoo  x  <  z  has  th«  constant 
wavs  nuMbsr  k(z  —  0}^  and  ws  dsnots  the  rsflsotsd  wavs  by 
u(z)  ik(2  —  0}(x*~z)j,  whsrs  ths  sxplloit  dspsndsnos 

of  ths  cosfflcisnt  of  rsflsotion^  ^i^)$  on  ths  position  of 
ths  left  faoe  of  the  slab,  z,  Is  indloatsd. 

Observe  that  when  z  »  0  ws  are  oonfrontsd  with  ths 
original  problem,  and  ^sn  z  »  b  ws  are  oonfrontsd  with  ths 
sliqple  problem  of  detemlning  ths  rsflsotion  oosffioisnt  for 
waves  incident  on  the  plans  interface  of  two  homogeneous 
media.  In  ths  latter  case  we  have 


(1) 


u(b) 


k(b  -  0)  -  kj^ 
k(b  -  0)  +  k^  ' 


which  follows  from  the  usual  assi^^ptions  that  are  made  con- 
csming  the  continuity  at  the  interface  of  ths  wave  function 
and  its  first  derivative  with  respect  to  z. 

We  now  wish  to  show  that  a  knowledge  of  the  reflection 
coefficient  for  the  Izihomogeneous  slab's  extending  from 
£  +  A  to  b  enables  one  to  find  ths  reflection  ooeffioient 
for  the  Blab  [z^b] .  This  is  aoeoKplished  through  use  of  the 
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localisation  principlo  given  in  [?] ,  valid  under  mild 
restriotiona  on  the  function  k(x)^  which  in  turn  was 
suggested  by  Bremmer's  paper,  jis]  «  Let  us  now  indicate  this 
principle . 

Return  for  a  moment  to  the  original  problem  with  z  »  0 
and  consider  a  plane  wave  expCik^x)  arriving  from  the 
homogeneous  space  7  <  0  and  entering  the  inhomogeneous 
slab.  At  the  botmdary  the  incident  wave  is  split  into  a 
reflected  wave  and  a  refracted  wave.  There  is  an  insoediate 
refliijoted  wave  and  an  immediate  refracted  wave  obtained  by 
supposing  that  the  inhomogeneous  medium  is  actually  homo¬ 
geneous  with  wave  nujnber  k(-{-  O).  If  the  inhomogeneous 
aediiuB  is  now  taken  to  be  the  limit  of  a  sequence  of  inter¬ 
faces  with  this  type  of  reflection  and  refraction  occurring 
at  each  interface,  we  can  obtain  the  total  reflected  wave 
and  the  disturbance  within  the  inhomogeneous  slab  by  adding 
up  the  effects  of  the  reflected  and  refracted  waves  obtained 
in  this  way  and  then  passing  to  the  limit.  This  principle  is 
valid  regardless  of  the  position  of  the  left  face  of  the  slab, 
z,  0  z  ^  b. 

To  apply  this  principle,  we  consider  the  situation 
illustrated  in  Fig.  3$  in  which  the  wave  numbers  of  the 
strata  are  as  shown. 
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Flg.  3 

Th«  O^osMtry  for  a  Class  of  Kodifisd  Problsns 

To  within  tonos  of  ordors  zsro  and  on«,  the  refleotion  oo* 
efficient  for  the  slab  j^z^bj  is  given  by  the  snm  of  three 
tenas.  The  first  is  the  isuediately  reflected  wave*  The 
second  is  the  wave  that  arlf  s  from  transaission  of  the 
incident  wave  through  the  interface  x  »  z,  refleotibn  at 
the  interface  x  +  d  and  transmission  through  the  inter- 
face  X  a  z.  Ihe  third  component  arises  from  trasismission  of 
the  incident  wave  throu^  the  interface  at  x  »  z,  reflection 
at  X  •>  z  +  d>  reflection  at  x  »>  z,  reflection  at 
X  «  z  +  d/  and  finally  transmission  through  the  interface 
at  X  »  z. 

MathematiosLlly,  we  are  led  to  the  equation 


u{z) 


k(g-Q)  ~  k:(g4-A-0) 
k(s-O)  +  k(z-fA-O) 
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+ _ __M5+fc5J _ 

k(sB-O)  +  k(2+A-0)  k(z-o)  +  k(z+A-0) 

k(a-O)  +  k(z+A-0) 

.  -  k(2r-0)^lk(z4-A~0)A^^^^^^^lk(z4-A~Q)A 

k(a-O)  +  k(z+A-0) 

- -  +  o(A). 

k(z+A-0)  +  k(a-O) 

A  passase  to  the  limit  in  which  A  tends  to  zero  then  shows 
that  u(z)  satisfies  the  Biccati  equation 

(3)  u'(*)  211cu-|^  u^, 

where  tbe  primes  denote  differentiation  with  respect  to  z. 

Integration  of  this  differential  equation,  which  is 
easily  aoooj^lished  numerically  on  a  hi^  speed  digital  com¬ 
puter  or  by  hand,  using  the  end  condition  of  equation  (1), 
on  the  interval  ^0,1^  then  yields  the  desired  reflection 
coefficient,  u(o). 

Similar  ELccati  equations  for  the  impedance  and 
reflection  coefficient  of  nonuniform  transmission  lines  have 
been  given  by  Pierce,  ^2lj  ,  and  by  Walker  and  Wax,  [26J  ;  see 
also  the  discussion  given  by  Schelkunoff ,  |^23j  ,  and 
Osterberg,  .  Ihe  idea  of  concentrating  attention  on  the 
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rftflftction  oooffioi«nt  as  a  funotioa  of  tho  thickness  of  tho 
slat)  goof  back  to  Ambarsonian,  [ij ,  who  applied  it  in  oertain 
radiative  transfer  problems  involving  diffuse  z*efleotion  from 
a  foggy  medium.  The  notion  was  oonsiderably  extended  and 
developed  by  Chandrasekhar j  [15]  and  Bellman  and  Kalaba,  [4]  . 

In  addition  it  forms  the  basin  for  the  treatment  by  Bellman, 
Kalaba,  and  Wing,  [8  -  i:^  of  various  problems  of  neutron 
multiplication  in  fissionable  material,  Inoluding  the  deter¬ 
mination  of  oritloal  mass. 

The  point  here  is  that  it  is  easier  to  deal  with  the 
nonlinear  lUooati  equation  and  the  end  condition,  insofar 
as  numerical  computation  is  concerned^  than  it  is  to  solve 
the  linear  two-point  boundary  value  problem  for  the  wave 
within  the  slab*  Furtherxsore,  In  numy  oases,  it  is  the 
reflected  wave  rather  than  the  wave  within  the  slab  which  is 
of  primary  physicaJL  ioportance.  Once  the  reflection 
coefficient  u(z)  has  been  determined,  though,  for  0  z  b, 
the  determination  of  the  wave  within  the  slab  is  also  reduced 
to  an  initial  value  problem. 

Similar  considerations  enable  us  to  determine  the  trans¬ 
mission  coefficient,  v(z).  For  z  »  b  the  problem  reduces 
to  that  of  determining  the  transmission  through  the  interface 
between  two  homogeneous  media  having  wave  numbers  k(b  —  0) 
and  kj^.  We  have  the  usual  end  condition 

(4)  v(b)  .  9)....  , 

k(b  -  0)  +  kj^ 
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To  oalculata  tht  transmission  ooeffiolont  for  the  case  in 
which  the  inho»ogen«o\i5  slab  extends  from  z  to  from  a 
knowledge  of  the  transmission  coefficient  for  the  case  in 
which  the  slab  extends  from  z  +  A  to  b,  we  write  the 
following  equation t 


(5)  r(z) 


2k(g-0) _ ^ik(  z+A-O )  A 

k{sM))  +  k(a+Ar-0) 


•»(z+4)  + 


•  k(z4-A-0)  -  k(a-O)^^ik(2H-A-0)A 
k(z+A-0)  +  k(2»-0) 


v(z+A)1  +  o(a). 


It  states  that  excluding  texms  of  order  higher  than  the  first 
in  A  the  total  transmitted  wave  consists  of  one  component 
idiich  is  transmitted  through  the  interface  x  »  z,  is  trans¬ 
mitted  through  the  Interface  x  «  z  -i-  A»  and  is  transmitted 
through  the  interface  x  «  b,  and  a  second  component  which 
arises  from  transsdssion  through  x  »  z^  reflection  at 
X  »  z  4-  A,  reflection  at  x  »  z  and  transmissions  at 
X  »  z  +  A  and  X  o  b.  A  passage  to  the  limit  in  which  A 
tends  toward  zero  then  ^elds  the  differential  equation 

(6)  "" 


AssuBiing  that  u(z}  has  been  calculated  on  the  interval 
0  ^  z  ^  b,  we  see  that  the  equations  (4)  and  (6)  permit  the 
determlDation  of  v(z)  on  the  same  interval j  5Ln  particular 
the  transmission  coefficient  for  the  inhomogeneous  slab,  v(0),, 
can  be  calculated. 

We  wish  to  point  out  the  strong  similarity  between 
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equatiozifi  (3)  and  (5)  ai>ov«  vhioh  desorlb*  rafleotion  and 
transnission  of  »av«s  and  the  eqnationa  in  (6.2)  of  our 
previoua  paper,  [s] ,  which  deacrlbs.  reflection  and  trana- 
mifaion  of  neutrona  in  a  rod  of  fisaionable  ssaterial. 
Abstractly  the  derivations  of  the  equations  for  the  wave  and 
the  paz*tiole  oases  are  identical.  We  plan  in  the  future  to 
investi^te  the  ioqplioations  of  this  wave-particle  analogy 
in  more  detail. 


4.  Stochastic  Wave  Propagation 

Let  us  now  turn  o\ir  attention  to  the  case  in  tdiioh  a 
wave  is  normal^  incident  on  a  stratified  slab,  each  stratum 
of  which  is  charaetezrlzed  by  a  wave  nuail^er  which  is  a  random 
variable.  Our  objectives  will  be  the  characterisation  of  the 
stochastic  reflected  and  refracted  waves. 


^0 

“a 

^1 

... 

|lefracted 

Wave 

h 

^2 

M 

"7 

Xq-O  Xg  X^g  X^j_-b 


gioident 

Wave 


Hsfleoted 
- ^ 

Wave 


Pig.  4 

A  Stochastic  Stratified  Slab  of  thickness  b 


Referring  to  the  above  figure,  we  denote  the  wave  number  of 
the  medium  z  <  0  by  and  that  of  the  stratum  between 
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Xgj  asid  by  the  randoa  Tarlable  a  =»  0,1,  2. 

For  B  >  »  b,  tht  medium  has  wave  aumbsp  kj^.  ®iu8  the 

slab  0  <  2  <  divided  into  N  —  1  strata.  We  put 

°  “  'Sh-I*  incident  wave,  which  arrives  at  the 

interface  z  »  from  +  oo,  is  specified  by  the 

function  e  ,  where  now  the  time  factor  e  is 

suppressed.  The  reflected  wave  is  denoted  by  R„e  , 

IV 

z  >  refracted  wave  is  denoted  by  Tj^e  , 


z  <  0. 


Our  attack  on  the  problem  of  the  determination  of  the 

random  variables  and  will  be  along  the  following 

lines.  We  first  imbed  the  problem  for  a  fixed  value  of  N 

within  the  class  of  problems  for  which  N  »  1,2,...  .  Next, 

using  the  localization  principle,  we  deid.ve  a  relationship 

between  the  random  variables  and  Lastly,  using 

analytic,  Monte  Carlo,  or  perturbation  techniques  we  determine 

successively  the  distribution  functions  for  the  reflection 

coefficient,  Rjj,  and  the  transmission  coefficient,  Tj^. 

In  this  paper  we  assume  that  the  random  wave  numbers  are 

independent,  but  no  a&9uxiq>tiona  of  smallness  of  variation  of 

the  random  nx2mber8  k^  about  their  means  need  be  made. 

The  sioqple  case  in  tdiich  H  *  1  is  easily  treated.  The 

ik,  E 

incoming  wave  is  e  ^  ,  z  >  0,  the  reflected  wave  is 
--lk,z  ik^z 

Rji^e  ^  ,  z  >  0,  and  the  refracted  wave  is  T^^e  ,  z  <  0. 

In  terms  of  the  wave  numbers  k^  and  k^,  the  reflection 

and  transmission  coefficients  sure  given  by 
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(1) 

(2) 


«  "  “o 

h-Tq-TT^^  • 


Th«8e  relations  follow  from  the  afts\isMd  continuity  of  tht  waTO 
funotion  and  its  first  dorivative  at  th«  intorfaco  z  ••  0. 

Wo  now  turn  to  the  case  illuatratod  in  Pig,  4.  total 
refleoted  wave  ic  oompostd  of  an  imaediatoly  ZMifleotod  wave^ 
a  wave  which  arises  from  transmisBion  through  the  intorfaoo 
at  z  m  reflootion  at  tho  face  z  *  3^2*  rttrani- 

mission  through  the  interface  at  z  »  ^ 

Expressed  mathematioally,  this  statement  becomes 


(3) 


idiere  the  dots  denote  the  contributions  from  terms  involving 
more  than  one  reflection  from  the  medium  z  <  To 

evaluate  these  terms  we  note  that  each  ooi^onent  incident  on 
the  interface  z  »  from  the  left  gives  rise  to  a  re¬ 

flected  wave  which  is  then  partially  reflected  at  the  inter¬ 
face  z  B  and  is  then  partially  transmitted  through  the 

interface  z  -«  ^Si^l*  coBqponents  are  the  tezms  of  a 

geometric  8eri«.a  with  coonon  ratio 


^1  ^  n  „’'^^1^N-1 

i^-n4 
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Equation  (3)  tooomos 


%  ^rncr  i  Z  7x  ,  ^ 


^  ,  ( ^Wl 


%-i 


%-l 


It  is  eonvoniont  to  introduce  the  local  reflection  and 
transffliasion  coefficients  for  waves  Incident  from  the  right 
on  the  interface  x  » 

/,.%  _  _  V-i 


and  the  corresponding  coefficients  for  waves  incident  fxvm 


the  left^ 


„  _  Vi 


4-Tc;~n^- 


Ihese  coefficients  obviously  satisfy  Stokes'  relations^ 


-  0, 


(10) 


-  1  -  r^. 


‘h  +  "  2. 


(11) 
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BO  that  aquation  (4)  becomes 

(«)  n,  .  r.  .  - :=®V?5 

Vi 


“H 


%-l 


r^+  e 


If«  finally,  we  make  use  of  Stokes'  relations  and  introduce 


(13) 


^H-1  “  • 


we  find 


(14) 


- -  ,  N«2,3,... 


+  r. 


hVi 


Equation  (1)  yields 


(15) 


h  “  'l' 


Equations  (14)  and  (15)  constitute  the  desired  stochastic 
functional  equations  for  the  reflection  ooeffloients 

Similar  considerations  for  the  transmission  coefficients 


yield  the  relations 


(16) 


V 


."‘^^l^K-l 


N 


N-2,3>..., 


Ti  .  t^. 


(17) 


W«  now  oonsi^isr  the  uses  to  whloh  we  rsar  put  the 
ftoohastlo  functional  equations  dezdved  in  the  last  section. 
Ihe  basic  problem  that  confronts  us  is  that  of  determining  the 
diBtzd.bution  of  and  given  the  distribution  of  the 

zrandoa  variables  r^  and  1  »  1,2, , .  .,N  —  1.  Here,  in 

speaking  of  the  detensination,  we  are  pri»arllj  Interested  in 
the  numerical  detemination. 

We  shall  discuss  three  possible  uses  of  the  preceding 
results . 

Monte  Carlo.  The  formula  seems  ideally  suited  for 
a  direct  determination  of  the  distribution  for  Rj^  and 
using  a  Monte  Carlo  technique  in  oonjunotion  with  a  digital 
ooB^uter  which  generates  the  x^ucdoa  sequences  ^r^|  and 
.  This  wlli  be  easy  to  carry  out  even  when  the  dlstri«- 
butlons  of  and  depend  xqpon  those  of 

^1—1*  runs  can  be  carried  out  in  a  natter  of 

minutes . 

2.  Reourreaoe  Relations  for  Dletributlon  Functions. 

Since  Rj^  is  ooBqples,  we  take  as  a  fundamental  distribution 
associated  with  this  stochastic  vaixlable,  the  function 

(l)  Pjj{s#r)  “  the  probability  that  ^  be  In  a  circle 

about  z  with  radius  r. 


Since  the  transformation  in  (4.14)  naps  the  interiors  of 
ci2rcles  into  the  interiors  of  circles,  it  is  clear  that  we 
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oan  obtain  a  rolation  for  pjj(z,r)  in  ttxnas  of  the  oorree- 
ponding  function  for  N  —  1,  [^t]* 

Since,  however,  p^(z»r)  depends  t^on  the  three  real 
variable*  x,  y  and  r,  whore  s  »  x  +  ly,  we  see  that 
this  relation  is  not  ideally  sxtited  for  direct  nusierioal 
cosiputation. 

Having  obtained  the  formula,  we  can  now  introduce 
perturbation  techniques*  Assuning  that  and  laay  be 
written  in  the  fom 

(2)  r^  «  r  + 

6^.  a* 

where  6  is  a  small  parameter  and  and  are  random 
variables  with  known  distributions,  we  can  writ* 

(3)  Pu(»»3?)  -  +  •••, 


and  use  the  recurrence  relations  to  obtain  relations  for 
Pjj  ^  ”h  f  •  •  •  • 

3.  Stochastic  Iteration,  The  formula  in  (4.l4)  yields 


a  relation  of  the  type 

(4) 


\  w^  +  ' 


¥ribitre  Ujj,  v^,  Wjj,  are  random  variables  determined  by  the 
sequences  |r^\  and 


As  a  aafc'bex*  of  fact,  know  that  (4.l4)  is  equivalent 


to  the  Bsatrix  relation 


(5) 


'"s 

V 

V  N 

1 

V  ’ 

J 

Vij 

These  relations  can  now  be  used  to  determine  the 
individual  and  Joint  aoogents  of  the  u^,  v,^, 
t*nm  of  the  noments  of  the  r^  and  6^,  of.  [2]  . 

Results  of  this  type  are  useful  in  oonneotion  with  the 
inverse  problem  of  detezioining  the  distribution  of  the  r^ 
and  given  the  obseirved  stOohastle  values  of  the 

Each  of  these  methods  will  have  their  analogues  in  the 
Kulti-dimenslonal  case.  From  the  standpoint  of  computational 
ease  and  effectiveness,  we  feel  that  the  Monte  Carlo  method 
is  the  Btost  promising. 


6.  Discussion 

It  is  clear  that  this  paper  is  but  a  first  step  toward 
a  more  comprehensive  theory  of  wave  propagation  in  stochastic 
media  based  on  the  principle  of  invariant  imbedding,  a 
radical  departure  from  current  practice.  In  future  papers,  we 
plan  to  discuss  applications  to  the  case  of  oblique  incidence, 
other  geometries  (see  [3/^o]  ),  and  electromagnetic  wave 
propagation,  per  se.  In  addition  we  propose  to  ejqploit  the 
wave-particle  similarity  mentioned  in  ^3* 

Finally,  let  us  mention  other  techniques  for  treating 
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nonlintair  equations  with  stoohastlo  eleiaents^  given  In 
Bellman,  [3] ,  iBuid  Kaiaba,  [l8] ,  and  whleh  are  relevant  In 
wave  propagation  studies. 
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